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Abstract 

A computer program is presented which calculates the elastic and inelastic scattering in inter- 
mediate and high energy nuclear collisions. A coupled-channels method is used for Coulomb and 
nuclear excitations of El, E2, E3, Ml, and M2 multipolarities, respectively. The program applies 
to an arbitrary nucleus, specified by the spins and energies of the levels and by reduced matrix ele- 
ments. For given bombarding conditions, the angular distribution of elastic and inelastic scattered 
particles and angular distributions of gamma-rays from the excited nucleus are computed. 

PACS numbers: 25.70.-z, 25.70.De 
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I. PROGRAM SUMMARY 



1. Title of program: DWEIKO (Distorted Wave EIKOnal Approximation) 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. 
Ireland 

Computers: The code has been created on an IBM-PC, but also runs on UNIX ma- 
chines. 

Operating systems: WINDOWS or UNIX 
Program language used: Fortran- 77 

Memory required to execute with typical data: 8 Mbytes of RAM memory and 1 MB 
of hard disk space 

No. of bits in a word: 64 

Memory required for test run with typical data: 1 MB 

No. of bytes in distributed program, including test data, etc. : 82238 

Distribution format: ASCII 

Keywords: Elastic scattering; Coulomb excitation; Relativistic collisions; Coupled- 
channels; Nuclear excitation 

Nature of physical problem: The program calculates elastic scattering differential cross 
sections, probabilities, and cross sections for inelastic scattering in nuclear collisions at 
intermediate and high energies {Eiab > 50 MeV/nucleon). It is particularly useful in 
the analysis of experiments with stable and unstable nuclear beams running at several 
intermediate-energy heavy ion accelerators around the world. 

Method of solution: Eikonal wavefunctions are used for the scattering. For each "im- 
pact parameter" entering the scattering matrix elements, one solves coupled-channels 
equations for the time-dependent Coulomb -|- nuclear field expanded into multipoles. 
A four-point Runge-Kutta procedure is used to solve the coupled-channels equa- 
tions. The elastic scattering is calculated purely with the eikonal approximation. The 
coupled-channels is a separate calculation for the inelastic amplitudes. The inelastic 
couplings, therefore, have no effect on the obtained elastic scattering cross sections. 
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Typical running time: Almost all the CPU time is consumed by the solution of the 
coupled-channels equations. It is about 2 min on a IGHz Intel P4-processor machine 
for the inclusion of 5 nuclear states. 

II. LONG WRITE-UP 
A. Introduction 

The eikonal approximation is very useful in the study of nucleus-nucleus scattering at high 
energies . In the Distorted Wave Born Approximation (DWBA) the transition amplitude 
for the reaction A{a, b)B involves a matrix element of the form p[ 

Tdwba = J ^^p^*{r){b,B\U,nt{r)\a,A)¥+\r)d'r^d\f,, (1) 

where Uint{i^) is the interaction potential, and \1/q, (^P/j) is the scattering wave function in 
the entrance (exit) channel, a = a + A{(3 = b + B). {a,A\ and (6, B\ are the initial and final 
intrinsic wavefunctions of the system, respectively. Using the eikonal approximation for the 
wave functions one has [|l| 

^(-)*(r)^(+)(r) ^exp{iq.r + zx(b)} , (2) 
where x(b) is the eikonal phase, given by 

1 

X(b) = -^J dz Uoptiv) . (3) 

The conditions of validity of the eikonal approximation are: (a) forward scattering, i.e., 
^ 1 radian, and (b) small energy transfers from the bombarding energy to the internal 
degrees of freedom of the projectile, or target. Both conditions apply perfectly well to direct 
processes in nuclear scattering at Eiah > 50 MeV/nucleon 

In the above equation, Uopt{'^) is the optical potential, with r = v^P~+^, where h can 
be interpreted as the impact parameter. For the Coulomb part of the optical potential this 
integral diverges. One solves this by using x = Xn + Xc, where xn is given by the equation 
above without the Coulomb potential and writing the Coulomb eikonal phase, xc as 

Xc{b) = 2r]\n{kb) , (4) 
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where r] = Z\Zie^ jlvf ^ Z\ and Z2 are the charges of projectile and target, respectively, v is 
their relative velocity, k their wavenumber in the center of mass system. Eq. ^ reproduces 
the exact Coulomb scattering amplitude when used in the calculation of the elastic scattering 
with the eikonal approximation [0: 

where 0o = ar 5fr(l +ir7/2). This is convenient for the numerical calculations since, as shown 
below, the elastic scattering amplitude can be written with the separated contribution of 
the Coulomb scattering amplitude. Then, the remaining integral (the second term on the 
r.h.s. of eq. ^4]) converges rapidly for the scattering at forward angles. 



Although the Coulomb phase in eq. ^ diverges at 6 = 0, this does not pose a real problem, 
since the strong absorption suppresses the scattering at small impact parameters. One can 
correct it for a finite charge distribution of the nucleus . For example, assuming a uniform 
charge distribution with radius R the Coulomb phase becomes 



Xc 



(b) = 2r] |e(6 - R) ln(fc6) + Q{R - h) \n{kR) + ln(l + ^/l - b^R^) 

-Vl- hVR' - \{l - h'/R'f"] } , (6) 

where G is the step function. This expression is finite for 6 = 0, contrary to eq. If one 
assumes a gaussian distribution of charge with radius i?, appropriate for light nuclei, the 
Coulomb phase becomes 

Xc{b)=2,^[\n{kb) + ]^E,{b'/R')] , (7) 
where the error function Ei is defined as 



00 



= I —dt. (8) 



X 



This phase also converges, as 6 ^ 0. The cost of using the expressions ^ and |^ is that 
the Coulomb scattering amplitude becomes more complicated than ^. Moreover, we have 
verified that the elastic and inelastic scattering cross sections change very little by using eqs. 
H or 1^, instead of eq. |[ 

The simplification introduced by the eikonal approximation is huge, as one avoids the 
calculation of scattering wavefunctions by solving numerically the Schrodinger equation for 
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each partial wave, as is done DWBA for nuclear scattering at low energies. For computer 
programs appropriate for nuclear scattering at low energies, see, e.g., the codes FRESCO 
§ and DWUCK4 §. 

The eikonal approximation is also valid in relativistic collisions, as it can be derived from 
relativistic wave equations, e.g., the Klein-Gordon equation [^]. Moreover, since it involves 
directions transverse to the beam, it is relativistically invariant. The eikonal wavefunction 
also allows the interpretation of the internal wave function variable b as an impact parameter. 
Thus, one can use the concept of classical trajectories, to obtain excitation amplitudes and 
use them as input to the scattering amplitudes. 

At intermediate energy collisions {Eiab — 50 MeV/nucleon), one must perform a correction 
due to the Coulomb deflection of the particle's trajectory. This correction amounts to 
calculating all elastic and inelastic integrals replacing the asymptotic impact parameter b 
by the distance of closest approach in Rutherford orbits, i.e., 

b' = ao + sjal + b^ , (9) 

where = ZiZ2e^/mv^ is half the distance of closest approach in a head-on collision of 
point charged particles. This correction leads to a considerable improvement of the eikonal 
amplitudes for the scattering of heavy systems in collisions at intermediate energies. 

As has been shown in ref. ||^, the nucleus-nucleus elastic and inelastic scattering at 
intermediate and high energies is corrected appropriately for relativistic kinematics if one 
replaces the quantity Oq in eq. |^ by Oq = ZiZ2e'^ /'jmv'^, where 7 = (1 — v^/c^)"-*^/^ is the 
Lorentz factor. 

III. THE OPTICAL POTENTIALS 

Usually, the optical potentials used in DWBA calculations are described in terms of 
Woods-Saxon functions, both for the real and for the imaginary part of the potentials, i.e., 

Uopt = -Vo fir, Rr, ar) - iW^ f{r, Ri, a^) , (10) 

where f{r,R,a) = 1/{1 + exp[(r — R)/a]}. The parameters entering these potentials are 
fitted to reproduce the elastic scattering data 0. 

In nuclear collisions at intermediate and high energies [Eiab > 50 MeV/nucleon) the 
elastic scattering data are very scarce. One has to resort to folding models with effective 
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interactions, at least as a guide for the experimental analysis. Among these models, the 
M3Y interaction is very popular. It has been shown to work quite reasonably for elastic and 
inelastic scattering of heavy ions at low and intermediate energy collisions |T0| . 

In its simplest form the M3Y interaction is given by two direct terms with different ranges, 
and an exchange term represented by a delta interaction: 

t{s) = A—- + B—- + C8{s), (11) 

PlS P2S 

where A=7999 MeV, B = -2134 MeV, C = -276 MeV fm^, A = 4 fm-\ and /^a = 2.5 
fm~^. The real part of the optical potential is obtained from a folding of this interaction 
with the ground state densities, pA and p^, of the nuclei A and B: 

UM3Y{r) = J dVi d^r2 PA(ri)pB(r2) t(s) , (12) 

with s = r + r2 — ri. The imaginary part of the optical potential is usually parameterized 
as Im Uopt = XUm3Y, with A = 0.6 - 0.8 g 0. 

This double folding M3Y potential yields values at the central region, r ~ 0, which are 
too large compared to usual optical potentials. However, one has to consider that the nuclear 
scattering at intermediate and high energies {Eiab > 50 MeV/nucleon) is mostly peripheral. 
Central collisions will lead to fragmentation reactions which are not being considered here. 
Thus, the only requirement here is that the optical potential reproduces well the peripheral 
processes. This is the case of the M3Y potential and the "t-pp" potential discussed below. 

Another simple method to relate the nuclear optical potential to the ground-state densi- 
ties is the "tpp" approximation. This approximation has been extensively discussed in the 
literature |jl2|, [T^. In its simplest version, neglecting the spin-orbit and surface terms, the 
optical potential for proton-nucleus collisions is given by 

Uoptir) =< tpn > Pn(r)+ < tpp > pp(r) (13) 

where p„ (pp) are the neutron (proton) ground state densities and < tpi > is the (isospin 
averaged) transition matrix element for nucleon-nucleon scattering at forward directions, 

^pi(q = 0) = fpi{q = 0) = — — api [ttpi + I) (14) 

where Opi is the free proton-nucleon cross section and is the ratio between the imaginary 
and the real part of the proton-nucleon scattering amplitude. The basic assumption here is 



that the scattering is given solely in terms of the forward proton-nucleon scattering amplitude 
and the local one-body density ||T^ . 

For nucleus-nucleus collisions, the extension of this method leads to an optical potential 
of the form 

Uoptij) = J < t^^(q = 0) > pA{r - r') psir') rfV , (15) 

where r is the distance between the center-of-mass of the nuclei. In this expression one uses 
the isospin average 



_ZiZ2 + NiN2. Z1N2 + Z2N1 

^ l^NN ^ — 1 1 '^DD H : : trm- 



(16) 



The parameters of the nucleon-nucleon cross scattering amplitudes for Eiab > 100 
MeV/nucleon are shown in Table II, extracted from ref. fl^. At lower energies one can 
use the isospin average values of Table II, which describes well the nucleus-nucleus elastic 



scattering at lower energies [TT 



E 


(^pp 


Oipp 


ipp 


^pn 


O-pn 


^pn 


[MeV/nucl] 


[mb] 




[fm2] 


[mb] 




[fm2] 


100 


33.2 


1.87 


0.66 


72.7 


1.00 


0.36 


150 


26.7 


1.53 


0.57 


50.2 


0.96 


0.58 


200 


23.6 


1.15 


0.56 


42.0 


0.71 


0.68 


325 


24.5 


0.45 


0.26 


36.1 


0.16 


0.36 


425 


27.4 


0.47 


0.21 


33.2 


0.25 


0.27 


550 


36.9 


0.32 


0.04 


35.5 


-0.24 


0.085 


650 


42.3 


0.16 


0.07 


37.7 


-0.35 


0.09 


800 


47.3 


0.06 


0.09 


37.9 


-0.20 


0.12 


1000 


47.2 


-0.09 


0.09 


39.2 


-0.46 


0.12 


2200 


44.7 


-0.17 


0.12 


42.0 


-0.50 


0.14 



Table I. Parameters ll^ for the nucleon-nucleon amplitude, as given by eq. |T7 



The formula |T5| can be improved to account for the scattering angle dependence of the 



nucleon-nucleon amplitudes. A good parametrization for the nucleon-nucleon scattering 
amplitude is given by 

kNN 



fNN{q) 



An 



-ctnn + (yNN) e 



(17) 



The nuclear scattering phase then becomes 



XN(h) 




dr dr' pi(r) 7ArAr(|b - s - s'|) p2{r') 



(18) 



where the profile function 7iVAr(b) is defined in terms of the two-dimensional Fourier trans- 
form of the elementary scattering amplitude 

7JV7v(b) = / exp [ - iq.h] fNN{q) dq , (19) 

2mkNN J L J 

and s, s' are the projections of the coordinate vectors r, r' of the nuclear densities on 
the plane perpendicular to the z-axis (beam-axis). For spherically symmetric ground-state 
densities eq. Ill reduces to the expression 



XN{h) 



dq q pA{q) fNN{q) Psiq) Jo{qh) , 



(20) 



where Pi{q) are the Fourier transforms of the ground state densities. 



Table II. Same as in tab 



E [MeV/nucl] 


< aNN > [fm^] 


< aNN > 


30 


19.6 


0.87 


38 


14.6 


0.89 


40 


13.5 


0.9 


49 


10.4 


0.94 


85 


6.1 


1 



e I, but for lower incident energies |11|. The values are averaged 



over pp and pn collisions. < > is taken as zero at these energies. 

The optical potential can also be obtained by using an inversion method of the eikonal 
phases. These phases might be chosen to fit the experimental data. In this approach, one 
uses the Abel transform 

hv d r x{b) 



Uoptir) 



mr dr (^52 _ j,2-)i/2 



r dr 



(21) 



This procedure has been tested in ref. leading to effective potentials which on the 

tail, where the process takes place, are very close to those obtained with phenomenological 
potentials of refs. |jl6|, |1^]. It can also be proved that under certain approximations, and for 



Gaussian density distributions, the potential obtained through eq. ^ coincides with that 
obtained with the double folding procedure |]15| . 
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IV. ELASTIC SCATTERING 



The elastic scattering in nucleus-nucleus collisions is a well established tool for the inves- 
tigation of ground state densities 0. This is because the optical potential can be related 
to the ground state densities by means of a folding of the nucleon-nucleon interaction with 
the nuclear densities of two colliding nuclei. But, as we have seen in the last section, this 
relationship is not straightforward. It depends on the effective interaction used, a proper 
treatment of polarization effects, and so on (for a review see, e.g., |T^]). At higher bombard- 
ing energies {E^ab > 50 MeV/nucleon), a direct relationship between the nuclear densities 
and the optical potential is possible, as long as the effects of multiple nucleon-nucleon scat- 
tering can be neglected [l^. The effects of real, or virtual, nuclear excitations should also 
be considered, especially for radioactive beams, involving small excitation energies. 

The calculation of elastic scattering amplitudes using eikonal wavef unctions, eq. H, is 
very simple. They are given by [ffl] 



Ue) = ik db b Uqb) |l - exp [zx(6)] } 



(22) 



where q = 2ksin{9/2), and 9 is the scattering angle. The elastic scattering cross section is 



dcTei 

'dn 



fern 



(23) 



For numerical purposes, it is convenient to make use of the analytical formula for the 
Coulomb scattering amplitude. Thus, if one adds and subtracts the Coulomb amplitude. 



fc{()) in eq. one gets 

fel{0) = fc{O)+tk 



db b Jo{qb) exp ixcip) 



exp 



Ml 



(24) 



where we replaced b in XNip) by b' as given by eq. ^ to account for the nuclear recoil, as 
explained at the end of section II. 



iXN{h) 



becomes zero for 



The advantage in using this formula is that the term 1 — exp 
impact parameters larger than the sum of the nuclear radii (grazing impact parameter). 
Thus, the integral needs to be performed only within a small range. In this formula, Xc is 
given by eq. | and fc{G) is given by eq. |^, with 



oo / 

j=0 ^-^ 



arctan ■ 



V 



J + 1 



(25) 



and C = 0.5772156... is the Euler's constant. 

At high energies the elastic scattering cross section for proton-nucleus collisions is also 
well described by means of the eikonal approximation [|1|] . The optical potential for proton- 
nucleus scattering is assumed to be of the form [R| 



where 



and 



Uoptir) = Uoir) + Us{r) (L.S) + Udr) 



Uo{r) = Vo fR{r) + iWo fi{r) + A i as Ws -^fi{r) 



— ) - :r.^-5 W 

m„c/ r dr 



(26) 



(27) 



(28) 



are the central and spin-orbit part of the potential, respectively, and Uc{r) is the proton- 
nucleus Coulomb potential. The Fermi (or Woods-Saxon) functions fi{r) are defined as 



before. The third term in eq. ^ accounts for an increase of probability for nucleon-nucleon 
collisions at the nuclear surface due to the Pauli principle. The existence of collective surface 
modes, the possibility of nucleon transfer between ions suffering peripheral collisions, and 
the possibility of breakup of the projectile and the target make possible enhanced absorption 
at the surface. The spin-orbit interaction in eq. ^ is usually parametrized in terms of the 
pion mass, 171^^: 2(^/m,rc)^ = 4 fm^. 

In the eikonal approximation, the proton-nucleus elastic scattering cross section is given 
by i 



daei 



F{e) + G{9) 



where 

F{9) = fc{9) + ik 
and 



dh h Jo{qb) exp ixcip) 



1 — exp 



ix{h) 



cos 



(29) 



khxs{h) (30) 



G{d) = ik dbb Ji{qh) exp ixc{b) + ^x(&) 



sm 



kh xs{b) 



(31) 



In the equation above q = 2ksm{9/2), where 9 is the scattering angle, x = Xn + Xc, -^o {Ji) 
is the zero (first) order Bessel function. The eikonal phase xs is given by |^ 



Xs{h) 



(32) 
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Eqs. p3|-^ describe the elastic scattering cross section of A (projectile) + B in the center 
of mass system. In the laboratory the scattering angle is given by |19[ 



where, p = Ma/ Mb, 



Ol = arctan <^ — — , > , (33) 

' 7[cos^ + p^(p,^i)] J 



Eiab [MeV/nucleon] 

C-l — ^ , [64:) 



where m-Ar is the nucleon mass, and 

7 is the relativistic Lorentz factor of the motion of the center of mass system with respect 
to the laboratory. 

The laboratory cross section is 

daei ,n N ^ {7^ [p^(p,gi) + cosg]^ + sin^g}^^^ '^(e\ (m^ 

V. TOTAL NUCLEAR REACTION CROSS SECTIONS 

The total nuclear inelastic cross section (including fragmentation processes) can be easily 
calculated within the optical limit of the Glauber model pi. It is given by 



an = 27 



poo 

In [1 - T{b)] b db , (37) 
Jo 

where T{b), the "transparency function", is given by 

Tib) = exp [2lmxN{b)] . (38) 

VI. THE SEMICLASSICAL METHOD AND COUPLED-CHANNELS PROBLEM 

Coulomb Excitation (CE) in high energy collisions is a well established tool to probe 



several aspects of nuclear structure |^ |21|, The CE induced by large-Z projectiles 

and/or targets, often yields large excitation probabilities in grazing collisions. This results 
from the large nuclear response to the acting electromagnetic fields. As a consequence, a 
strong coupling between the excited states is expected. 
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Since there will be very little deflection by the Coulomb field in collisions with impact 
parameter greater than the grazing one, the excitation amplitudes can be calculated assum- 
ing a straight-line trajectory for the projectile. A small Coulomb deflection correction can 
be used at the end, with the recipe given by eq. ^. 

We describe next a method for the calculation of multiple excitation among a finite num- 
ber of nuclear states. The system of coupled differential equations for the time- dependent 
amplitudes of the eigenstates of the free nucleus is solved numerically for electric dipole 
(El), electric quadrupole (E2), electric octupole (E3), magnetic dipole (Ml), and magnetic 
quadrupole (M2) excitations. 

Similarly, one can also calculate the amplitudes for (nuclear) monopole, dipole, 
quadrupole, and octupole excitations. 

In high energy nuclear collisions, the wavelength associated to the projectile-target rela- 
tive motion is much smaller than the characteristic lengths of the system. It is, therefore, a 
reasonable approximation to treat r as a classical variable r(t), given at each instant by the 
trajectory followed by the relative motion. At high energies it is also a good approximation 
to replace this trajectory by a straight line. The intrinsic dynamics can then be handled 
as a quantum mechanics problem with a time-dependent Hamiltonian. This treatment is 
discussed in full details by Alder and Winther in ref. ||23|| . 

The intrinsic state {ipit) > satisfies the Schrodinger equation 

{Ho + V [r{t)]} I m) = ^^^-^^ • (39) 

Above, Hq is the intrinsic Hamiltonian and V is the channel-coupling interaction. 

Expanding the wave function in the set {| j); j = 1,A^} of eigenstates of Hq, where N 
is the number states included in the coupled-channels (CC) problem, we obtain a set of 
coupled equations. Taking the scalar product with each of the states < k\, we get 

N 

ih dk{t) = ^ {k\ V{t) I j) exp [i{Ek - Ej)t/h] a^it) , A; = 1 to , (40) 

where is the energy of the state \n) . It should be remarked that the amplitudes depend 
also on the impact parameter h specifying the classical trajectory followed by the system. 
For the sake of keeping the notation simple, we do not indicate this dependence explicitly. 
We write, therefore, a„(t) instead of a„(6, t), restoring the notation with 6, or t, whenever 
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necessary. Since the interaction V vanishes as t ^ ±00, the amphtudes have as initial 
condition a„(t —00) = 6ni and they tend to constant values as t ^ 00. 

A convenient measure of time is given by the dimensionless quantity r = '~fvt/b, where 
7 = (1 — v^/c^)^^/^ is the Lorentz factor for the projectile velocity v. A convenient measure 
of energy is Eq = 'yhv/b. In terms of these quantities the CC equations become 



dak{T) .v^ /, I X I .\ /-^ N / N Ti// \ ^(^) Ek-Ej 
W{t) I j) exp {likjT) aj-(r) ; W{t) = ; - 



(41) 

The nuclear states are specified by the spin quantum numbers I and M. Therefore, 
the excitation probability of an intrinsic state | n) =| In,Mn) in a collision with impact 
parameter b is obtained from an average over the initial orientation (Mi) , and a sum over 
the final orientation of the nucleus, respectively: 

^nW = ^ E \<]mM'- (42) 
The total cross section for excitation of the state |n > is obtained by the classical expression 

a„ = 27r / Pn{b) b db . (43) 



VII. TIME-DEPENDENT ELECTROMAGNETIC INTERACTION 

We consider a nucleus 2 which is at rest and a projectile nucleus 1 which moves along 
the z-axis. Nucleus 2 is excited from the initial state \IjMj > to the state \IkMk > by the 
electromagnetic field of nucleus 1. The nuclear states are specified by the spin quantum 
numbers Ij, Ik and by the corresponding magnetic quantum numbers Mj and M^. We 
assume that nucleus 1 moves along a straight-line trajectory with impact parameter b, which 
is therefore also the distance of the closest approach between the center of mass of the two 
nuclei at the time t = 0. The interaction, Vc{t), due to the electromagnetic field of the 
nucleus 1 acting on the charges and currents nucleus 2 can be expanded into multipoles, as 



explained in ref. PH|. One has 



Wc{r) = y^ = J2w.x,{T) , (44) 



ttA/j. 
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where tt = E, M denotes electric and magnetic interactions, respectively, and 



1 Zie 1 



2n 



' ' hvb^ A y (2A + 
where A^(7rA,/i) is the multipole moment of order A/i P^j, 



(45) 



X(EA,/i) = J dh pc(r) Fi^(r) 



(46) 



and 



A^(Ml,/i) = -^ / c/V Jc(r).Lfryi, 



i.J , (47) 

Pc (Jc) being the nuclear charge (current). The quantities QT,\^{r) were calculated in ref. 
P3| , and for the El, E2, and Ml multipolarities. We will extend the formalism to M2 
and E3 excitations. 



We use here the notation of Edmonds |23] where the reduced multipole matrix element 
is defined by 



(48) 



To simplify the expression (^) we introduce the dimensionless parameter V'lj'^'' by the 
relation 

,A+i Zie 1 



^g'^^ = (-1) 



hvh^ A 



y^^^f^-M..(vrA,-/.). 



Then we may write eq. ^ in the form 

N 



dokir) 
dr 



r=l ttA/^ 



(49) 



The explicit expressions for eq. ^ can also be obtained by a Fourier transform of the 
excitation amplitudes found in ref. [^, i.e., 

W^^,{t) = Tr-7r V^x,iuj)duj = -- / e^«- KAM(0^e, (50) 

where u = {Ek — Ej)/h = Eq^/H (here we omit the sub- indexes kj for convenience). The 
expressions for K-A/x('-^ji) are given by ||2^ 



KrAfi(^^ij) 



V7 



-1)'^v/2ATT(^) G.A^(^)ir^(6,)-M(vrA,-/i). (51) 
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Using the properties of KrAf((t^) for negative uj^ one can show that KrA^(— 

Tit I 



1^— 1)^+^V^*^ (cij). Then, one gets from eq. EO 



^"'^^^^ = 2^ ^"^^'^ (jeT ^^2ATTG.am (^) FA^(6„r)>((7rA,-/i), (52) 



27r/iv7 

where 



/•oo 

(e, r) = 2 / cos (^r) (0 rf^, for A + /i = even, 
Jo 

= - / sin (er) e'i^^ (0 c^e, for A + /x = odd. (53) 
^ Jo 

These integrals can be obtained analytically. 



Using the functions Gj^x^ (c/v) derived in ref. we get explicit closed forms for the El, 
E2, E3, Ml and M2 multipolarities, 

QElo{^^J, r) = -^r(l)\T) ■ QEl±l{^^J, t) = T^^t) , (54) 

QMioi^ij, r) = 0; QAn±ii^ij, t) = i (^^ j 0'^(r) , (55) 

QE2oiC^„r) = VQ{2r^-l)<P\r) , 
QE2±l{^^J, t) = ±3 (2 - P^) 7r0^(r) ; Qe2±2{t) = 30^(r) , (56) 

QM2o{C,ij, r) = , 

QM2±l{^^J, r) = -3i (^) 7r0^(r) ; QM2±2{r) = tSi (^) 0^(r) , (57) 



and 



QE3o{^^J, r) = -3^ ^7 (5 - f3') r {2t' - 3) 0' (r) , 
QE3±li^^J, t) = (15 - 11(3') 7^ {At' - l) 0^(r) 



QE3±2{^^,, r) = -7y (3 - P') r0^(r); gi,3±3(6„ r) = T^^^'i^) 



(58) 



where (r) = (1 + r^)"^''^, and (3 = v/c. 

The fields QnXfii^ij, t) peak around r = 0, and decrease rapidly within an interval Ar ^ 1, 
corresponding to a collision time At 6/7V. This means that numerically one needs to 
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integrate the CC equations in time within an interval of range n x At around r = 0, with 
n equal to a small integer number. 

It is important to notice that the multipole interactions derived in this section assume 
that the Coulomb potential has an 1/r shape (when 7 — > 1) even inside the nuclei. This 
is a simplification justified by the strong absorption at small impact parameters for which 
the Coulomb potential deviates from its 1/r form. We have not included the deviations 
to the 1/r behavior of the Coulomb interaction inside the nuclei, as the proper relativistic 
treatment of it is rather complicated. This has been discussed in details in ref. 2^ . 



VIII. TIME-DEPENDENT NUCLEAR EXCITATION: COLLECTIVE MODEL 

In peripheral collisions the nuclear interaction between the ions can also induce excita- 
tions. According to the collective, or Bohr-Mottelson, particle-vibrator coupling model the 
matrix element for the transition j — > k is given by H, [2^ 



4M,|V)v(A;.)|/,M, >= -^7|= < 4M,|nj/,M, > ^^(r) U,{r) (59) 

where 6\ is the vibrational amplitude and Ux{r) is the transition potential. To follow the 
convention of ref. we use ao instead of 6q in the equation above. 

The transition potentials for nuclear excitations can be related to the optical potential 



in the elastic channel. This is discussed in details in ref. p9[. The transition potentials for 
isoscalar excitations are 

Uo{r) = 3UoAr) + r^^^ , (60) 

for monopole. 



, , , ^ dU,„,t . 1 ^ d'U, 
Ui{r) = 

for dipole, and 



ar 

for quadrupole and octupole modes. Rq is the nuclear radius at | the central nuclear density. 

The deformation length 6x can be directly related to the reduced matrix elements for 
electromagnetic transitions. Using well-known sum-rules for these matrix elements one finds 
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a relation between the deformation length, and the nuclear sizes and the excitation energies. 



For isoscalar excitations one obtains 

1 



a, 







27r 



A (2A + 1) -— 

6 rriN ALx 



(63) 



where A is the atomic number, < > is the r.m.s. radius of the nucleus, and is the 
excitation energy. 



For dipole isovector excitations 29 



71 A 1 



(64) 



2 niN NZ E^ ' 

where Z [N) the charge (neutron) number. The transition potential in this case is modified 
from eq. ^ to account for the isospin dependence |^ . It is given by 

'N - Z^ 



UAr) 



-A 



A 



dr 3 ° dr"^ 



(65) 



where the factor A depends on the difference between the proton and the neutron matter 
radii as 



A 



2{N -Z) 
3l 



R-n Rp 
2 {.Rn ~l~ Rp 



AR 



■np 



Rn 



(66) 



Thus, the strength of isovector excitations increases with the difference between the neutron 
and the proton matter radii. This difference is accentuated for neutron-rich nuclei and 
should be a good test for the quantity ARnp which enters the above equations. 

Notice that the reduced transition probability for electromagnetic transitions is defined 
by H 



B(7rX:i 



3] 



1 



2/. + 1 

These can be related to the deformation parameters by 



</,-||A^(7rA)||/, >| 



5(E0) 



212 



and 



3Zei?i 
lOvr 



B{EX) 



a. 



' 



B{E1) 



9 fNZe\ 



An 



A>2 



Air 



-ZeR 



A-l 



(67) 



(68) 



(69) 



The time dependence of the matrix elements above can be obtained by making a Lorentz 
boost. One gets 



-7 



V2A + 1 



< hM,\Yx,\IjMj > Y,^i9it),0) Ux[rit)] 



(70) 
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where r{t) = ^JW^^^)W^ = hj (^{r\ 9 = T(j){T), and 



< IkMk\Yx^\I,M, > 



-r 



(24 + l)(2A + l)l'/' 



Jfc A Ij 
-Mk Ai Mj 



Ik A Ij 





• (71) 



To put it in the same notation as in eq 
coupled-channels equations become 



4vr(2J, + 1) 

1), we define Q);i;(r) = V'^'l^^^it) / E^, and the 



da 



N 



dr 



/ J / J / J I ^N\^i 

j = l X/l TT 



^ 52 52 52 I ^) + Qc^Xf^i^kj, t) tpif exp (i^kjT) ajir) . (72) 



IX. ABSORPTION AT SMALL IMPACT PARAMETERS 



If the optical potential ?7opt(r) is known, the absorption probability in grazing collisions 
can be calculated in the eikonal approximation as 



A(b) = exp 



' 2 f°° 

— J lm[Uoptir)]dz 



(73) 



where r = y/b"^ + z^. If the optical potential is not known, the absorption probability can 
be calculated from the optical limit of the Glauber theory of multiple scattering (also from 
the "t-pp" approximation), which yields: 



A{b) = exp i^-aNN j j pi(r') P2(r - r') rfV dz^ . 



(74) 



where cttvat is the nucleon-nucleon cross section and pi is the ground state density of the 
nucleus i. For stable nuclei, these densities are taken from the droplet model densities of 
Myers and Swiatecki but can be easily replaced by more realistic densities. 

Including absorption, the total cross section for excitation of the state |n > is obtained 

by 

(Tn = 2-K j A{h)Pn{h) hdh . (75) 
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X. ANGULAR DISTRIBUTION OF INELASTICALLY SCATTERED PARTI- 
CLES 



The angular distribution of the inelastically scattered particles can be obtained from the 
semiclassical amplitudes, af^^j^^{b), described in section VII. For the excitation of a generic 
state |n >, it is given by 



/•oo 

fLiiO) = ik dbb Uqb) e^^(^) a,ib) , (76) 

^0 

where we simplified the notation: = aY^]Mn^ with /i = M„ — Mi. 

The inelastic scattering cross section is obtained by an average over the initial spin and 
a sum over the final spin: 

do'inel 



^ E \fLi\' ■ (77) 



The program DWEIKO uses eq. ^ to calculate the angular distribution in inelastic scat- 
tering. But it is instructive to show how it relates to the usual semiclassical approximation. 
For collisions at high energies, the integrand of eq. oscillates wildly at the relevant impact 
parameters and scattering angles. One can use the approximation 



together with the stationary-phase approximation 

n / r, • \ 1/2 

j G(x)e''^(")dx ~ [riij^] G'(a;o)e''^("°) , (79) 
where Xq is the point of stationary phase, satisfying 

0'(xo) = . (80) 
This approximation is valid for a slowly varying function G{x). 



Only the second term in the brackets of eq. |7^ will have a positive {b = bo > 0) stationary 



point, and eq. becomes 
fl^,{e) ^ I— v^exp [Imx^(6o)] exp Mbo) + m(m + l/2)/2] a^{bo) , (81) 
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where 

^=-qb + 2ri \n{kb) + Rbx'n (b) , (82) 
and bo, the "classical impact parameter" is the solution of 

-g + -^ + Rex'^(&o) = 0. (83) 

This equation has 2 solutions: (a) one corresponding to close (or nearside) collisions, (b) 
and another corresponding to far (or farside) collisions. These are collisions passing by one 
side and the opposite side of the target, but leading to the same scattering angle. They thus 
lead to interferences in the cross sections. 

In collisions at high energies, the inelastic scattering is dominated by close collisions and, 
moreover, one can neglect the third term in eq. PBI. The condition 0'(&o) = implies 

We observe that the relational is the same [with cot(^/2) ~ sin^"'^(^^/2)] as that between 
the impact parameter and the deflection angle of a particle following a classical Rutherford 
trajectory. 



Inserting these results in eqs. 76 and O, one gets 



dn 



(^)27TT ^ l<:M„(^o)re--(M. (85) 
\ q / ZJi + i ^^^^^^^ 

One can easily see that the factor Arfk'^/q'^ is the Rutherford cross section. 

The above results show that the description of the inelastic scattering in terms of the 
eikonal approximation reproduces the expected result, i.e., that the excitation cross sec- 
tions are determined by the product of the Rutherford cross sections and the excitation 
probabilities. This is a commonly used procedure in Coulomb excitation at low energies. 

The cross sections in the laboratory system are obtained according to the same prescrip- 
tion as described at the end of section IV. 



XI. ANGULAR DISTRIBUTION OF 7-RAYS 

After the excitation, the nuclear state |//) can decay by gamma emission to another state 
\Ig)- Complications arise from the fact that the nuclear levels are not only populated by 
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(a) 



(b) 



FIG. 1: Schematic description of a nuclear excitation (solid line) followed by 7-decay (solid wavy 
line), (a) The dashed lines are transitions due to internal conversion (unobserved). The dashed 
wavy line is an unobserved 7-decay. (b) Direct emission of an observed gamma ray. 

Coulomb excitation, but also by conversion and 7-transitions cascading down from higher 
states (see figure |l](a)). To compute the angular distributions one must know the parameters 
A; {i — > i) and e/ {i — > j) for / > 1 [^ , 



(i — , j) =ai{i — > j) {i — > j) 



(86) 



where ai is the total /-pole conversion coefficient, and 



A, 



ttI 



Stt (/ + 1) 



21+1 



1/2 



{2Ij + l)-^/^ {lj\\f^'^Minl)\\k) 



(87) 



_l[{2l + l)\\f h Vc. 

with s{l) = I for electric (tt = E) and s{l) = Z + 1 for magnetic (vr = M) transitions. The 
square of A^^i is the /-pole 7-transition rate (in sec~^). 

As for the non-relativistic case P3|, 153], the angular distributions of gamma rays following 
the excitation depend on the frame of reference used. In our notation, the z-axis corresponds 
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to the beam axis, and the statistical tensors are given by (we use the notation of ||23|, ^ 



(88) 

Ml ^ 

where / is the state from which the gamma ray is emitted, and 1 denotes the initial state 
of the nucleus, before the excitation. To calculate the angular distributions of the gamma 
rays one needs the statistical tensors for A; = 0, 2, 4 and —k < n <k (see [p3| , p3|). 

Instead of the diagram of figure 0(a), we will consider here the much simpler situation in 
which the 7-ray is emitted directly from the final excited state / to a lower state which 
is observed experimentally (see figure |l|(b)). The probability amplitude for this process is 

ai^f^g = J2 ihMgka \H,\ IjMf) , (89) 

where {IgMgha \ H^ \ IfMf) is the matrix element for the transition f g due to the emis- 
sion of a photon with momentum k and polarization a. The operator accounts for 
this transition. The angular dependence of the 7-rays is given explicitly by the spherical 
coordinates 6 and of the vector k. 

Since the angular emission probability will be normalized to unity, we can drop constant 
factors and write it as (an average over initial spins is included) 



E«W {I,Mgka\H,\IfMf) 



Mf 



(90) 



The transition operator can be written as 

ff, = E^r = E^Jr^®A^:^^ (91) 

Im Im 

where the first operator in the sum acts between nuclear states, whereas the second operator 
acts between photon states of well defined angular momentum, l,m. 

Expanding the photon state |kcr) in a complete set |Zm) of the photon angular momentum, 
and using the Wigner-Eckart theorem (angular momentum notation of ref. P^), one gets 



l.m 



-1)^/--^EI ^ \{ka\lm){lg\mif). (92) 
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One can rewrite Ikcx) in terms of |zcr), i.e., in terms of a photon propagating in the z- 
direction. This is accomphshed by rotating |ko") to the z-axis, using of the rotation matrix 

(ka|/m) = J^D'mm' k) {za\lm') . (93) 

m' 

Expanding the photon field in terms of angular momentum eigenfunctions, one can show 
that m, M 



'^d^m for Ti = E 
(zal/mvr) = <( Vl ^ (94) 

^(t5,„ for 7r = M. 

One has now to express the operator O^,^, in eq. |9l| in terms of the electric and magnetic 
multipole parts of the photon field. This problem is tedious but straightforward [^ ). 
Inserting eqs. ^ and ^ in eq. yields (neglecting constant factors) 



(J,M,ka \H,\ IjMf) = J2 i-lY^-'^''^' J {21 + 1) {21 J + 1) 



Lm 



Mf m ~Mf 



X D^^^ (z k) [A^, + a Ami] , (95) 



where A^/ is given by eq. 87 



Inserting eq. ^ into eq. |90| one gets a series of sums over the intermediate values of the 
angular momenta 

/ If 
Mg m -Mf 



M^,Mg,cr,Mj, 



^(_^)M,.MH-' M^ ] DljDiJ A,A^, , (96) 

m' -M'f J L J 

where A; = A^;^ + IAmi- The product A;A*, is always real since (— l)**^'^ = 11 (the parity). 

Assuming that the particles are detected symmetrically around the z-axis one can inte- 
grate over (pparticie, what is equivalent to integrating, or averaging, over 0^. This yields the 
following integral 

V v47r \ m -m \ a -a 



23 



To simplify further eq. ^ we use (see ref. ||23|, p. 441, eq. II. A. 61) 




Mg m' -M'j 



/ V k 



m —m' K 



If k 



I V 
^ If 



and 



I J I' 



I J I' 



<T=(-1,1) 



a -a 



A^iAli, = I \ 1 -1 

0, 

V 

where use has been made of the parity selection rule 



A^/A;^, , for j = even 
for j = odd, 



1) , for electric transitions 



-1) , for magnetic transitions . 
Eq. ^ becomes 

Wie)= Yl (-1)'™'+'+'"^ ai^fal^f,sj{2l + 1) {21' + 1) {21 j + 1) 

1,1' ,m,m' 

' If If k\ (if I' k 

Mf -M'f K j yMf -m' K 



X (2j + 1) {2k + 1) 



X 



I J V \ I J V 



I I' k 



10-1/ \ m -m' j \ I'f If Ig 



A/A;,Pj(cos^) . 



Using 



I J v \ln I' 

m -m' / \ 1 -1 



l+l'+k 



1 



2k + 1 



one gets 



E (-iy'-'^V^2/ + l)(2/' + l)(2/^ + l)(2fc + l)|a._;r 



A; — even, Af^ , 

Mf,l,l' 



If If k\ I J V 



I I' k 



Mf -Mf ) \1 -1 I \ I'f If Ig 



A;A;*Pj (cos( 
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or in a more compact form 



1,1' 



(100) 

where 

Fk (/, r, //) = (-1)'^"'^-' ^(2/+l)(2/' + l)(2/; + l)(2A;+l) 

/' k\ j I r k] 
[i -I o) [If If i,j ' 

The angular distribution of 7-rays described above is in the reference frame of the excited 
nucleus. To obtain the distribution in the laboratory one has to perform the transformation 

f sin^ 1 

Ol = arctan <^ — } , (102) 



, 7 [cos 9 + 13] 
and 

W{eL) =1^1 + (3 cos efW (6) , (103) 
where 7 is given by eq. and /3 = ^Jl — 1/7^. The photon energy in the laboratory is 



A. Computer program and user's manual 

All nuclear quantities, either known from experiments or calculated from a model, as well 
as the conditions realized in the experiment, are explicitly specified as input parameters. 
The program DWEIKO then computes the optical potentials (if required), differential cross 
section for elastic scattering, and Coulomb + nuclear excitation probabilities and cross 
sections, as well as the angular distribution of the 7-rays. 

The units used in the program are fm (femtometer) for distances and MeV for energies. 
The output cross sections are given in millibarns. 



1. Input parameters 



To avoid exceeding use of computer's memory, the file DWEIKO. DIM contains the dimension of 
the arrays and sets in the maximum number of levels (NMAX), maximum total number of magnetic 

25 



substates, (NSTMAX), maximum number of impact parameters (NBMAX), and maximum number 
of coordinates points used in the optical potentials and absorption factors, (NGRID). A good 
estimate is NSTMAX = {2Jmax + 1)NST, where Jmax is the maximum angular momentum of the 
input states. 

Most integrals are performed by the 1/3-Simpson's integration rule. It is required that NGRID 
be a even number, since an extra point (origin) is generated in the program. 
The input file allows for comment lines. These should start with a '7^' sign. 
The file DWEIKO.IN contains all other input parameters. These are 

1. AP, ZP, AT, ZT, which are the projectile and the target mass and charge numbers, respec- 
tively. 

2. EGA, the bombarding energy per nucleon in MeV. 

3. EX(j) and SPIN(j): the energy and spins of the individual states j. 

4. MATEl(j,k), MATE2(j,k), MATES (j,k), MATMl(j,k), MATM2(j,k) the reduced matrix 
elements for El, E2, E3 and Ml, M2 excitations, j — > (as defined in (^8|)), in e fm 
(El, Ml), e fm^ (E2,M2), and e fm'^ (E3) units. 

5. F(0,j), F(l,j), F(2,j), F(3,j), the fractions of sum rule of the deformation parameters for 
monopole, dipole, quadrupole, and octupole nuclear excitations, entering eq. (|70|). 

To simplify the input, the deformation parameters are calculated internally in DWEIKO 
using the sum rules 6^56, with Ex replaced by the energy of the corresponding state, EX(j). 



The user needs to enter the fraction of those sum rules exhausted by the state j, i.e., the 
program uses 5'^ = f\5\, with (0 < /a < 1) entered by the user, and 5\ given by eqs. 63-3^. 



The input cards in file DWEIKO.IN are organized as following: 

1. AP, ZP, AT, ZT, EGA 

Gharges and masses (AP,ZP,AT,ZT), bombarding energy per nucleon in MeV/nucleon. 

2. IW, lOPM, lOELAS, lOINEL, lOGAM 
IW=0(1) for projectile (target) excitation. 
IOPM=1(0) for output (none) of optical model potentials. 
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IOELAS=(0)[1]2 for (no output) [center of mass] laboratory elastic scattering cross section. 

IOINEL=(0)[1]2 for (no output) [center of mass] laboratory inelastic scattering cross section. 

IOGAM=(0)[1]2 for (no output) [output of statistical tensors] output of gamma-ray angular 
distributions. 

The statistical tensors are calculated for each impact parameter, so that one can use them in 
the computation of dP^N^Mib) /d^l.y = PM{h).dW^M^M / dTt^. If IW = 0, a transformation 
to the laboratory system is performed. 

3. NB, ACCUR, BMIN, lOB 

NB = number of impact parameter points (NB < NBMAX). 

ACCUR = accuracy required for the time integration of the CC-equations for each impact 
parameter. A reasonable value is ACCUR = 0.001, i.e., 0.1%. 

BMIN=minimum impact parameter (enter for default. The program will integrate from 
6 = 0, with strong absorption included). 

IOB=1(0) prints (does not print) out impact parameter probabilities. 

4. lOPW, lOPNUC 

lOPW is a switch for the optical potential model (0PM). 

IOPW=0 (no OMP, IOELAS=0), 1 (Woods-Saxon), 2 (read), 3 (t-pp folding potential), 4 
(M3Y folding potential). 

IOPNUC=0 (no nuclear), 1 (vibrational excitations). 

If the optical potential is provided (I0PW=2), it should be stored in 'optw.in' in rows of 
R X Real[U(R)] x Imag[U(R)]. The program makes an interpolation to obtain intermediate 
values. 

The first line in 'optw.in' is the number of rows (maximum=NGRID) . 

5. VO [MeV], rO [fm], d [fm], VI [MeV], rOJ [fm], dl [fm] 

If I0PW=1, enter VO [VI] = real [imaginary] part (>0) of Woods-Saxon potential, rO [rOJ] 
= radius parameter (R = rO * (APV^ + AT^/^), d [dl] = diffuseness. 

If lOPW is not equal to 1, place a sign at the beginning of this line, or delete it. 
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6. VSO [MeV], rO_S [fm], dS [fm], V_surf [fm], a_surf [fm] 

If I0PW=1 and AP, or AT, equal to one (proton), enter here spin-orbit part. If not, place 
a sign at the beginning of this line, or delete it. 

VSO = depth parameter of the spin-orbit potential (>0), rO_S = radius parameter, dS = 
diffuseness, V_surf = depth parameter of the surface potential (>0), a_surf = diffuseness, 

7. Wrat 

If I0PW=4, enter Wrat = ratio of imaginary to real part of M3Y interaction. 

If lOPW is not equal to 4, place a sign at the beginning of this line, or delete it. 

8. THMAX, NTHETA 

If I0ELAS=1,2 or I0INEL=1,2 enter here THMAX, maximum angle (in degrees and in the 
center of mass), and NTHETA (< NGRID), the number of points in the scattering angle. 

If JOEL AST or lOINEL are not 1, or 2, place a sign at the beginning of this line, or 
delete it. 

9. JINEL 

If I0INEL=1 enter the state (JINEL) for the inelastic angular distribution (1 < JINEL < 
NST) . 

If lOINEL is not 1, or 2 place a sign at the beginning of this line, or delete it. 

10. NST 

NST (< NMAX) is the number of nuclear levels. 

11. I, EX(I), SPIN(I) 

Input of state labels (I), energy EX(I), and angular momentum SPIN (I). I ranges from 1 to 
NST and should be listed in increasing value of energies. 

12. I, J, El[e fm], E2[e fm^], E3[e fm^], Ml[e fm], M2[e fm^] 

Reduced matrix elements for El, E2, E3, Ml and M2 excitations: < Ij\\0{E /M\L)\\Ii >, j 
> i , for the electromagnetic transitions. Matrices for reorientation effects, i— >^i, can also be 
given. 
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Add a row of zeros at the end of this hst. If no electromagnetic excitation is wanted just 
enter a row of zeros. 

13. J, FO, Fl, F2, F3 

If I0PNUC=1 enter sum rule fraction of nuclear deformation parameters for monopole, 
dipole, quadrupole nuclear excitations (ALPHA0,DELTE1,DELTE2,DELTE3) for each ex- 
cited state J. 

If IOPNUC=0 insert a comment card ('#') in front of each entry row, or delete them. 

14. IFF, IGG, THMIN, THMAX, NTHETA 

If I0GAM=2, enter IFF,IGG = initial and final states (IFF > IGG) for the gamma transi- 
tion. 

THMIN, THMAX = minimum and maximum values of gamma-ray angles (in degrees) in 
the laboratory frame. 

NTHETA = number of angle points (< NGRID). 
2. Computer program 

The program starts with a catalogue of the nuclear levels by doing a correspondence of 
integers to each magnetic substate. J = 1 corresponds to the lowest energy level, with the 
magnetic quantum number Mi = —I^. J increases with Mi and so on for the subsequent 
levels. 

A mesh in impact parameter is done, reserving half of the impact parameter points, 
i.e., NB/2, to a finer mesh around the grazing impact parameter, defined as ho = 
fm. The interval bo/ 3 fm < b < 2bo fm is covered by this mesh. A 
second mesh, with the other half of points, extends from b — 2bo fm to b — 200 fm. Except 
for the very low excitation energies {Ej <C 1 MeV), combined with very large bombarding 
energies (7 3> 1), this upper value of b corresponds to very small excitation probabilities, 
and the calculation can be safely stopped. The reason for a finer mesh at small impact 
parameters is to get a good accuracy at the region where both nuclear. Coulomb, and ab- 
sorption factors play equally important roles. At large impact parameters the probabilities 
fall off smoothly with b, justifying a wider integration step. 
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A mesh of NGRID points in polar coordinates is implemented to calculate the nuclear 
excitation potentials and absorption factors, according to the equations presented in sections 
2.2 and 2.3. The first and second derivatives of the optical potentials are calculated by the 
routine DERIVATIVE. A 6-point formula is used for the purpose. The routines TWOFOLD 
computes the folding over the densities, as used in eq. (|73p. Routines RHOPP and RHONP 
generate the liquid drop densities, and the routine PHNUC computes the eikonal integral 
appearing in eq. (|73|). 

Repeated factors for the nuclear and for the Coulomb potentials are calculated in the 
main program and stored in the main program with the arrays PSITOT and PSINUC. 
These arrays are carried over in a common block to the routine VINT which computes the 
function Q^^i^i^kj, t) + QcnX/ii^kj, r), used in eq. (^. 

The optical potentials are read in routine OMP_READ, or calculated in routines 
OMP_WS (Woods-Saxon), OMPJDEN (t-pp) and 0MP_M3Y (M3Y). Routine TM3Y sets 
the M3Y interaction. DEFORM calculates the effective potentials in the Bohr-Mottelson 
model. SIGNNE and PHNNE return the nucleon-nucleon cross sections and the parameters 
of the nucleon-nucleon scattering amplitude PHNUCF calculates the eikonal phase in the 
"t-pp" approximation with the help of the Fourier transform of the ground state densities, 
provided by FOURIER. 

The time integrals are performed by means of an adaptive Runge-Kutta method. All 
routines used for this purpose have been taken from the Numerical Recipes [^. They 
are composed by the routines ODEINT, RKQS, RKCK, and RK4. The routine ODEINT 
varies the time step sizes to achieve the desired accuracy, controlled by the input parameter 
ACCUR. The right side of ( [72D is computed in the routine DCADT, used externally by the 
fourth-order Runge-Kutta routine RK4. RKCK is a driver to increase time steps in RK4, 
and RKQS is used in ODEINT for the variation of step size and accuracy control. The main 
program returns a warning if the summed errors for all magnetic substates is larger than 10 
X ACCUR. 

Elastic scattering is calculated within the routine ELAST for nucleus-nucleus, and 
ELASTP for proton-nucleus, collisions. Inelastic scattering is calculated within the rou- 
tine INELAST, and 7-ray angular distributions are calculated in the routine GAMDIS. 

The routine THREE J computes Wigner-3J coefficients (and Clebsh-Gordan coefficients), 
RACAH the 6j-symbols, or Racah coefficients, YLM is used to compute the spherical har- 
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monies, LEGENGC the Legendre polynomials, and BESS JO (BESSJl) [BESSJN] the Bessel 
function Jq (-^i) [Jn]- 

The routines SPLINE and SPLINT perform a spline interpolation of the excitation am- 
plitudes, before they are used for integration by means of the routines QTRAP and QSIMP, 
from Numerical Recipes PB[ . 



XII. TEST INPUT AND THINGS TO DO 

A test input is shown below. It corresponds to the excitation of giant resonance states 
in Pb by means of the reaction ^^O (84 MeV/nucleon) + ^''^Pb. Assuming that an isolated 
state is excited, and that it exhausts fully the sum rules, one gets {B{E\) = B{E\, E^)) 

B(El) = A JL ^ e^ (104) 

and, for (A > 1) 

3 n \ Z'^ I A, for isoscalar excitations: 

B{E\) = ^ A(2A + 1) e^x { ' ' (105) 



47r E^ 2m 



N 



NZ/A, for isovector excitations; 



From these values, the reduced matrix elements can be calculated from the definition in eq. 

m 

Example 1 - The input list below calculates the excitation cross sections for the isovector 
giant dipole (IVGDR) and isoscalar giant quadrupole (ISGQR) resonances in Pb, for the 
above mentioned reaction. The user should first run this sample case and verify that the 
output numbers check those in ref. [^. In particular, compare the results for elastic 



scattering with the upper figure 2 of ref. [^. Also compare the inelastic excitation of the 



IVGDR with the upper figure 3 of ref. . Following this it might be instructive to change 
the input of energies, spins, and excitation strengths for low lying states, giant resonances, 
etc. 

# Input of program 'DWEIKO' 

# Ap Zp At Zt Einc[MeV/u] 
17 8 208 82 84. 

# IW=0(1) IOPM=0(1) IOELAS=0(1)[2] IOINEL=0(1)[2] IOGAM=0(1)[2] 
11112 
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# NB ACCUR BMIN[fm] IOB=1(0) 
200 0.001 0. 

# lOPW lOPNUC 
1 1 

# VO [MeV] rO[fm] d[fm] VI [MeV] rOJ [fm] dl [fm] 
50. 1.067 0.8 58. 1.067 0.8 

# VSO [MeV] rO_S [fm] dS [fm] Vsurf [fm] dsurf 

# 15. 1.02 0.6 50. 0.8 

# Wrat 

# 1. 

# THMAX NTHETA 
6. 150 

# JINEL 
3 

# NST 
3 

# I Ex[MeV] SPIN 

1 

2 10.9 2 

3 13.5 1 

# i -> j El[e fm] E2[e fm**2] E3[e fm**3] Ml[e fm] M2[e fm** 2] 
1 2 76.8 

1 3 9.3 


# j FO Fl F2 F3 (NOTE: fractions of sum rules for deformation parameters) 

2 0.6 

3 1.1 

# IFF IGG THMIN THMAX NTHETA 
2 1 20. 70. 150 

Example 2 - Table III gives the results of an experiment on Coulomb excitation of S and 
Ar isotopes Choose the minimum impact parameter, BMIN, so that it reproduces the 
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maximum scattering angle (^max = 4.1°) in the experiment (use the formula h = cot(^/2) 
and read the last paragraph of section II). Using the B{E2) values in the table reproduce 
the cross section values (careful with the units!). 



Secondary beam 




40g 


42g 


44 


46 


Eiab [MeV/nucleon] 


39.2 


39.5 


40.6 


33.5 


35.2 


Energy of the first excited state [MeV] 


1.286(19) 


0.891(13 


0.890(15) 


1.144(17) 


1.554(26) 


cr(E2;0+,, -.2+; ^lab > 4.1°) [mb] 


59(7) 


94(9) 


128(19) 


81(9) 


53(10) 


BiE2; 0+,. - 2+) [e^ fm^] 


235(30) 


334(36) 


397(63) 


345(41) 


196(39) 



Table III. Experimental results on Coulomb excitation of S and Ar projectiles impinging 

on a ^^^Au target [^. 



Example 3 - Compare the results of this code with those from the ECIS code ^ ^ . 
But notice that the ECIS code has relativistic corrections in kinematic variables only. The 
relativistic dynamics in the Coulomb and nuclear interaction are not accounted for. Thus, 
one should expect disagreements for high energy collisions {Eiab > 100 MeV/nucleon). 



XIII. OUTPUT 

The output of DWEIKO are in the files 

1. DWEIKO. OUT: Probabilities and cross sections; 

2. DWEIKO.OMP.OUT: Optical model potential; 

3. DWEIKO_ELAS.OUT: Elastic scattering cross section; 

4. DWEIKO-INEL.OUT: Inelastic scattering cross section; 

5. DWEIKO_STAT.OUT: Statistical tensors; 

6. DWEIKO_GAM.OUT: Angular distributions of gamma-rays. 
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